transZeta = function(zeta){
	res = zeta
	if (zeta <= 1){
		z = qchisq(zeta, 2)
		print(z)
		res = 0.006*z + 0.0005
	}
	res = res/10
}

x <- matrix(seq(0, 1, 0.001), ncol=1)
y <- apply(x, MARGIN=1, transZeta)
plot(x, y, type="l", main="Ponderation de zeta", col="red", las=1, xlab="zeta", ylab="f(zeta)")

